(B) AMENDMENTS TO THE CLAIMS 



1. (Original): A method for seismic data processing, comprising: 

constructing explicit depth extrapolation operators with variable operator lengths 
depending on maximum dip angle, accuracy condition, and wavenumber; and 
constructing operator tables using the explicit depth extrapolation operators. 

2. (Original): The method of claim 1, further comprising: 

performing depth migration using the explicit depth extrapolation operators from the 
operator tables. 

3. (Original): The method of claim 1, wherein the step of constructing operator tables 
comprises: 

selecting a maximum dip angle; 
selecting a type of operator; 

selecting an accuracy condition for the type of operator; 
selecting a plurality of wavenumbers; and 

performing the following steps for each of the plurality of wavenumbers: 
selecting a plurality of operator lengths for the selected wavenumber; and 
performing the following steps for each of the plurality of operator lengths: 

determining if the selected type of operator with the selected operator length 
satisfies the selected accuracy condition at the selected maximum dip angle at the 
selected wavenumber; 

determining if the operator length is the smallest operator length satisfying the 
accuracy condition at the selected maximum dip angle for the selected 
wavenumber; and 

storing the operator length in an operator table if the operator length is the 
smallest operator length satisfying the accuracy condition at the maximum dip 
angle for the selected wavenumber. 
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4. (Currently amended): The method of claim 3, wherein the steps performed for each 

wavenumber cbmprise: 

selecting a plurality of pairs of first operator half-lengths and second operator half- 
lengths for the selected wavenumber; and 

performing the following steps for each of the plurality of pairs of first operator half- 
lengths and second operator half-lengths: 

determining if the selected type of operator with the selected pair of first operator 
half-length and second operator half-length satisfies the selected accuracy condition 
at the selected maximum dip angle at the selected wavenumber; 
determining if the s e l e ct e d typ e of op e rator with th e first operator half-length is the 
smallest first operator half-length satisfying the accuracy condition for the selected 
wavenumber and the second operator half-length is the smallest second operator half- 
length satisfying the accuracy condition for the selected wavenumber; and 
storing the pair of first operator half-length and second operator half-length in an 
operator table if the first operator half-length is the smallest first operator half-length 
satisfying the accuracy condition for the selected wavenumber and the second 
operator half-length is the smallest second operator half-length satisfying the 
accuracy condition at the selected maximum dip angle for the selected wavenumber. 



5. (Currently amended): The method of claim 1, wherein the step of constructing operator 
tables comprises: 

selecting a maximum dip angle; 

selecting a type of operator; 

selecting an accuracy condition for the type of operator; 

determining preliminary operator tables for varying with variable operator lengths; 
selecting a plurality of wavenumbers; and 

performing the following steps for each of the plurality of wavenumbers: 

selecting a smallest operator length in the preliminary operator tables satisfying the 
accuracy condition for the selected wavenumber; and 
storing the smallest operator length in an operator table. 
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6. (Currently amended): The method of claim [[4]] 5, wherein the steps performed for each 
wavenumber comprise: 

selecting a pair of smallest first operator half-length and smallest second operator half- 
length in the preliminary operator tables satisfying the accuracy condition for the selected 
wavenumber; and 

storing the pair of smallest first operator half-length and smallest second operator half- 
length in an operator table. 

7. (Currently amended): The method of claim 2, wherein the step of performing depth 
migration comprises: 

selecting a seismic data set; 

selecting a plurality of depths in the seismic data set; 
selecting sub-domains for each of the plurality of depths; 
determining a velocity model for the seismic data set and the plurality of depths; 
determining a lowest velocity for each sub-domain from the velocity model; 
selecting an operator table with varying variable operator lengths; 
selecting a plurality of frequencies in the selected sub-domain; and 
performing the following steps for each of the plurality of frequencies: 
performing the following steps for each of the plurality of depths: 

performing the following steps for each of the sub-domains at the depth: 

calculating a highest wavenumber in the sub-domain from the frequency and 

the lowest velocity for the sub-domain; 

selecting a maximum operator length for the calculated highest wavenumber 
from the operator table; 

applying operators from the operator table with operator lengths no greater 
than the maximum operator length to the selected sub-domain. 

8. (Currently amended): The method of claim 7, wherein the operator table with varying 
variable operator lengths comprises explicit extrapolation operators with varying variable 
operator lengths. 
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9. (Currently amended): The method of claim 7, wherein the step of selecting an operator table 
with varying variable operator lengths further comprises: 

interpolating the operator table with varying variable operator lengths. 

10. (Original): The method of claim 7, wherein each operator length comprises a pair of first 
operator half-length and second operator half-length. 

11. (Original): The method of claim 3, wherein the step of selecting a type of operator comprises 
applying the following formula: 

P(x i9 y j9 a> 9 z + Az) = 

Z X >ym>K (*/ > yj > z )> p ( x i - x i > yj -y m >*>> z l 

where P(x if y j} co, z+Az) is a seismic wavefield at lateral location jc z = / Ax, y } = /Ay, and depth 
z+Az; co is angular frequency; P(x it y Jf co, z) is a seismic wavefield in a sub-domain £) at depth z; 
Ax, Ay, and Az are step lengths in the x-, y- and z-coordinate directions, respectively; 
LfkJ^fD, z)) and M(k m max (D, z)) are operator half-lengths in the x- and y- coordinate directions, 
respectively, for the highest wavenumber ka, max (D, z) in the sub-domain D at depth z; W(x u yj, 
kcofc, y jt z), Az) is an explicit depth extrapolation operator; and k m (x u yp z) is a local wavenumber. 

12. (Original): The method of claim 11, wherein the wavenumber k m (xu yj, z) is given by the 
formula: 

K( x i>yj> z ) = — e r. 

c(x i9 y J9 z) 
where c(x u y jt z) is a local propagation velocity. 
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13. (Original): The method of claim 11, wherein the highest wavenumber k^fD, z) for a 
frequency co for a sub-domain D at depth z is given by the formula: 

C(Az)= " , 

where c min (D, z) is the lowest velocity for the sub-domain D at depth z. 

14. (Original): The method of claim 3, wherein the step of step of selecting a type of operator 
comprises applying the following formula: 

P(x i , <y, z + Az) = w ( x i > K (*« > z )> p ( x i -x l9 a> 9 z), 

where P(x u co, z+Az) is a seismic wavefield at lateral location x t = i-Ax and depth z+Az, co is 
angular frequency, P(x i} co, z) is a seismic wavefield in a sub-domain D at depth z, Ax and Az are 
step lengths in the x- and z-coordinate directions, respectively, L(kco max (D, z)) is an operator half- 
length in the x-coordinate direction for the highest wavenumber kco max (D, z) in the sub-domain D 
at depth z; W(x u k m (x u z), Az) is an explicit depth extrapolation operator, and k w (x it z) is a local 
wavenumber. 

15. (Original): The method of claim 14, wherein the wavenumber kjx u z) is given by the 
formula: 

c(x n z) 

where c(x if z) is a local propagation velocity. 
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16. (Original): The method of claim 3, wherein the step of selecting a type of operator comprises 
applying the following formula: 

P(x ( 9 y j9 co 9 z + Az) = Y< W i (*. > yj > z )> A*) h i 0/ > Vj 

/=o 

where P(x it y Jf co, z+Az) is a seismic wavefield at lateral location x t = i'Ax 9 yj = j Ay, and depth 
z+Az; co is angular frequency; P(x iy y jt co, z) is a seismic wavefield in a sub-domain D at depth z; 
Ax, Ay, and Az are step lengths in the x-, y- and z-coordinate directions, respectively; 
LflcJ^fD, z)) is an operator length for the highest wavenumber kco max (D, z) in the sub-domain D 
at depth z; Wi(k m (xu y jt z), Az) is an explicit depth extrapolation operator expressed in cylindrical 
coordinates; k m (x if y jf z) is a local wavenumber; and hi(x it y jt co, z) are auxiliary fields resulting 
from applying recursive 2D Chebyshev filters to the seismic wavefield P(x u yj, co, z). 
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